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CD ■ Abstract 

We study equilibrium and non-equilibrium aspects of the normal state of cold 
and dense QCD and QED. The exchange of dynamically screened magnetic 
^ \ gluons (photons) leads to infrared singularities in the fermion propagator for 

excitations near the Fermi surface and the breakdown of the Fermi liquid de- 
scription. We implement a resummation of these infrared divergences via the 
Q^ . Euclidean renormalization group to obtain the spectral density, dispersion 

f-^ \ relation, widths and wave function renormalization for single quasiparticles 

^D \ near the Fermi surface. We find that all feature scaling with anomalous di- 

^ \ mensions: u)p{k) ex |A; - A;i?| i^^2a ; r(A;) oc |A; - kF\'^^^^ ; Zp{k) (x\k - kp\'i^^^ 

I ; with A = ^ for QED , ^^^ for QCD with N^ colors and Np flavors. 

Q \ The discontinuity of the quasiparticle distribution at the Fermi surface van- 

^. : ishes. For k ^ kp ^e find n^^^^ = ^ig^ _ ^^^^^ + 0{k - kp? with M 

the dynamical screening scale of magnetic gluons (photons). The dynamical 
^ . renormalization group is implemented to study non-equilibrium relaxation. 

The amplitude of single quasiparticle states with momentum near the Fermi 
surface falls off as |^fc^fc^(t)| ^ |*fc^fc^(to)|e"'^('')(*"*«)[to/t]^^- Thus quasi- 
particle states with Fermi momentum have zero group velocity and relax with 
a power law with a coupling dependent anomalous dimension. 

I. INTRODUCTION AND MOTIVATION 

There is a substantial theoretical and experimental effort to map the phase diagram of 
QCD as a function of temperature (T) and chemical potential (/i). Current theoretical ideas 
suggest |jl|] that heavy ion collision experiments from SIS to AGS, SPS, RHIC and finally 
LHC have the potential to study the region of the phase diagram for T < 300 — 400 Mev 
and /i < 0.6 Gev with higher (T) and lower (yu) for RHIC and LHC. Understanding this 



^ 



region of the phase diagram can provide insights into the QCD phase transition in the 
Early Universe, about 10 /is after the Big Bang as well as the equation of state of hot 
and dense QCD. Matter at low temperature < 10 Mev and up to nuclear matter density 
Po ~ 0.16 fm~^ ; fj, ^ 300 Mev is amenable of study by low energy nuclear systems such 
as multifragmentation phenomena in nuclei. Cold-dense nuclear matter for densities larger 
than a few times nuclear matter density cannot be studied with terrestrial accelerators and 
is the realm of astrophysical compact objects, such as neutron stars [@,|1. The fascinanting 
possibility of detecting a phase transition in quark matter in neutron star X-ray binaries 
was raised recently [||,H, where the signal would be a pronounced peak in the frequency 
distribution of X-ray neutron stars due to a long spin-up stage and the cooling history 
as revealed by the (soft) X-ray spectra p2|- p^] . Thus while QCD at high temperature 
and relatively small chemical potential can be experimentally studied with ultrarelativistic 
heavy ion collisions, astrophysical observation of the properties of neutron stars can provide 
observable signatures from cold and dense QCD if quark matter is the correct description of 
the core of spinning neutron stars. 

While there is a substantial body of results on QCD at finite temperature on the lattice, 
the lack of manifest reality of the fermion determinant with finite chemical potential presents 
a problem for the lattice program with the important exception of two colors which has 
received considerable attention recently ^-^. 

An important aspect of QCD at large chemical potential is that of color superconductivity 
01 which arises from a paring instability of the free Fermi gas in the color antitriplet channel. 
Since the original proposal of color superconductivity via a one (screened) gluon exchange 
0, there has been an increased interest in color superconductivity |T0[- ^^ and diquark 
condensation [^]. The presence of diquark condensates in the cold and dense core of neutron 



stars could have potential observable signatures in their cooling history as well as in the 
magnetic fields of pulsars p2|- p^] . Therefore the study of cold-dense QCD is warranted by a 
definite phenomenological and observable impact if not on terrestrial accelerator experiments 



certainly in the astrophysical signatures of neutron stars p5[| . 

Goals: The common framework to study degenerate correlated fermion systems is that 
of Fermi liquid theory (see next section). The emergence of superconductivity (diquark 
condensation in the case of QCD) is associated with the instability of the normal Fermi 
liquid towards an attractive pairing interaction. In the case of a weakly interacting Fermi 
system the starting point is the free Fermi gas and pairing results in the opening of a gap in 
the single particle spectrum at the Fermi surface. Fermi liquid theory is argued to describe 



the low energy effective field theory of nuclear matter p6[ and is therefore an important tool 



to study nuclear matter and its impact in astrophysical compact objects. 

Recently it has been argued, within the context of color superconductivity, that the 
exchange of dynamically screened (via Landau damping) magnetic gluons results in strong 
infrared divergences that lead to the breakdown of the Fermi liquid description of cold-dense 
QCD in perturbation theory |^. A similar situation was found in the case of a non- 



relativistic electron gas with magnetic interactions pB|-p(][|, where it was argued that there 
would be no observable consequences of the breakdown of Fermi liquid theory for terrestrial 



densities. 

While diquark condensation and color superconductivity in its various forms have been 
studied extensively in the literature, we are aware of only one previous study of non-Fermi 
liquid aspects of cold dense QCD [^, which revealed large corrections to the superconducting 
gap |2^. Further corrections to the superconducting gap were found from lifetime effects 
||3T| not associated with the breakdown of Fermi liquid theory. 

Studying Fermi liquid aspects of the normal phase is an important part of the program 
towards understanding cold-dense QCD. As mentioned above, understanding the properties 
of the normal phase is perhaps the first step towards a complete assessment of the pairing 
instabilities and properties of the superconducting state. Furthermore, if the pairing in- 
teraction opens a gap at the Fermi surface of some quarks, such as the two or three color 
superconducting phases (2SC or 3SC) [|TU|- ]TB| , P(]| , pB| , the remaining gapless quarks will be 
described by the concomitant Fermi liquid. 

Our goal in this article is to provide a comprehensive study of non-Fermi liquid aspects 
in the normal phase of cold-dense QCD postponing to a forthcoming article the study of the 
implications of the breakdown of Fermi liquid theory on color superconductivity. While the 
study in |^ focused on the corrections to the color superconducting gap and issues of gauge 



invariance, we study both equilibrium and non-equilibrium aspects of the non-Fermi liquid 
behavior. In particular we focus on i) dispersion relation and damping rates of quasiparticles 
near the Fermi surface, these reveal anomalous dimensions resulting from the breakdown of 
Fermi liquid theory and ii) the relaxation of these quasiparticles: again we find anomalous 
relaxation with origin in the same infrared divergences responsible for the breakdown of 
Fermi liquid theory. Our main motivations for initiating this study and long term goals are 
manifold: a) to obtain a further assessment of non-Fermi liquid corrections to the super- 
conducting gaps, critical temperature and spectrum of excitations in the superconducting 
phase, b) a study of the potential implications of non-Fermi liquid corrections to the neu- 
trino emissivity and the cooling rate of neutron stars with degenerate quark-matter cores, 
c) a more complete and detailed understanding of the properties of dense QCD in a regime 
which is not yet amenable to lattice simulations, d) a study of transport phenomena in the 
non- Fermi liquid, which is relevant to cooling and thermodynamics of neutron stars. 

In this article we begin this program by studying in detail the breakdown of the Fermi 
liquid description and by providing a comprehensive analysis of the spectrum of single quasi- 
particle excitations in the normal phase along with their relaxation properties. 

Strategy: We study both QED and QCD at zero temperature but large (baryon) den- 
sity so that a perturbative analysis is reliable. In this regime the hard-dense-loop (HDL) 
approximation, which is the finite density equivalent of the hard-thermal-loop program of 
Braaten and Pisarski for finite temperature [p2|-|37[], is reliable and describes the main as- 
pects of static and dynamical screening of gluons (and photons) . The leading order in the 
HDL approximation is the same in the abelian (QED) and the non-abelian (QCD) theories 
and screening of gauge fields is completely determined by the one loop quark polarization 
at finite (and large) density pi|-p7[]. In this approximation, static "longitudinal" gluons (in- 
stantaneous Coulomb interaction) are screened by a Debye screening mass mu oc gfi while 



transverse gluons are only dynamically screened via Landau damping ||3^J3^ . To this order in 
the HDL approximation the polarization tensor for gluons in QCD is similar to that for pho- 
tons in QED save for trivial color and flavor factors. The main difference between QCD and 
QED in this approximation is that while one gluon exchange leads to an attractive (pairing) 
interaction in the antitriplet particle-particle channel and therefore to diquark condensation, 
there is no such attractive channel in QED0. Thus to this order in the HDL approximation, 
the Fermi liquid aspects of the normal state of cold-dense QCD are similar to those of QED. 
Thus we present our study within the framework of QED which accounting for the trivial 
color and flavor factors also describe those of the normal state of QCD to leading order in 
the HDL approximation. 

As described in detail in section II, a Fermi liquid description has an associated "order 
parameter" , this is the jump discontinuity of the Fermi distribution function (of the inter- 
acting) system at the Fermi momentum. This discontinuity is given by the residue (wave 
function renormalization) of the quasiparticle pole for quasiparticles with the Fermi momen- 
tum. The breakdown of the Fermi liquid picture is associated with the vanishing of this 
order parameter, i.e, the Fermi distribution function is continuous at the Fermi momentum. 
We begin by obtaining the quark propagators to leading order in the HDL approximation, 
corresponding to one bare gluon exchange in the quark self-energy and show explicitly that 
to this order there is a sharp discontinuity at the Fermi surface determined by the wave 
function renormalization of the quasiparticle pole. 

However this picture does not survive screening corrections to the gluon propagator, 
whereas longitudinal gluon exchange leads to an infrared finite contribution, the exchange 
of magnetic gluons which are only dynamically screened by Landau damping introduce log- 
arithmic divergences in the quark propagator for quasiparticles near the Fermi surface. 

For excitations with Fermi momentum we argue that these infrared divergences are sim- 
ilar to those of a critical theory at the upper critical dimensionality. Thus we provide a 
resummation of the quark propagator for particles with the Fermi momentum using the Eu- 
clidean renormalization group which reveals the emergence of anomalous dimensions in the 
spectral density. Non-equilibrium aspects are studied by implementing a dynamical renor- 
malization group IQ which provides a resummation of the quark propagator directly in real 



time. The dynamical renormalization group reveals power law relaxation with anomalous 
dimension for quasiparticles with Fermi momentum. 

Summary of the results: The exchange of dynamically screened magnetic gluons 
leads to infrared divergences in the single particle propagator for particle excitations near 
the Fermi surface. These divergences are akin to those found in critical phenomena for a 
critical theory at the upper critical dimensionality. We implement a resummation of the 
perturbative expansion via the Euclidean renormalization group. We find that the particle 



^However we expect the Overhauser effect to be present in QED |19], i.e, particles above and 
holes below the Fermi surface bound by their mutual (screened) Coulomb attraction. These are 
the counterpart of exciton bound states in condensed matter. 



component of the quark propagator for excitations near the Fermi surface is a scaling function 
of the two variables u = u — fi ; k = k — /i with anomalous exponents that depend on the 
gauge coupling. For k ^ the spectral density has the Breit-Wigner (quasiparticle) form, 
with the quasiparticle dispersion relation, width and residue given respectively by 



ujp{k) oc |A;|i-2^ ; T{k) oc sin[7rA]|A;| 1-2^ ; Zp{k) oc |A;|i-2^ 
with the effective coupling given by 
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for QCD with A''^ colors and Np flavors. 



'1.21 



The residue of the quasiparticle pole vanishes as A; — > fc^, leading to the following form 
of the (quasi) particle distribution function near the Fermi surface 



nkRikf 



sin[7rA] 
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f for QED 

f ^ for QCD with Np flavors. 
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revealing the vanishing of the discontinuity of the Fermi distribution function at the Fermi 
surface and therefore the vanishing of the Fermi-liquid order parameter. 

Implementing a real-time, i.e, a dynamical version of the renormalization group to study 
the non-equilibrium relaxation of single quasiparticles near the Fermi surface. We find that 
the amplitude of the wave function of single quasiparticle states near the Fermi surface fall 
off as 



\^k^kAt)\ oc |^^...,(to)| e-^(^)(*-*°) 
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Thus quasiparticles with the Fermi momentum have vanishing group velocity and relax 
with an anomalous power law. 

We make a comparison between these features of cold-dense QCD (and QED) and those 
of strongly correlated quasi-one dimensional Fermi systems that feature non-Fermi liquid 
behavior and are described as Luttinger liquids p9| , ^ . 



The article is organized as follows: in section II we summarize the aspects of Fermi 
liquid that are relevant for our discussion. In section III we obtain the expression for the 
quasiparticle distribution function in terms of the quark spectral density and obtain the 
equation of motion for quarks which will be used to study non-equilibrium aspects. In section 
IV we study the equilibrium aspects of single quasiparticles. We begin by studying the quark 
propagator to lowest order in the HDL approximation, i.e, with the self energy given by the 
exchange of hard (bare) gluons and make contact with the Fermi liquid description to this 
order. Soft {q < gfi) gluons require HDL resummation, and the propagator for particles near 
the Fermi surface is computed by including HDL (screening) corrections to the exchanged 
gluon. The resulting infrared divergences are recognized to be similar to those of a critical 
theory at its upper critical dimension and resummed using the Euclidean renormalization 



group. The renormalization group improved spectral density features scaling behavior that 
lead to the single quasiparticle dispersion relation and lifetime that scales with anomalous 
dimensions. We show that the jump discontinuity of the single particle distribution function 
vanishes at the Fermi surface as a consequence of the vanishing of the single quasiparticle 
residue at the Fermi momentum. Section V explores non-equilibrium aspects: the relaxation 
of single quasiparticle excitations near the Fermi surface. Implementing a real-time version of 
the renormalization group reveals that single quasiparticle states with Fermi momentum relax 
with a power law with anomalous dimension. In section VI we summarize the connection 
between QED and the normal state of QCD to the order studied and address the important 
issue of vertex corrections. Our conclusions are discussed in section VII. In this section we 
also discuss the striking resemblance of the spectral density and relaxation to that obtained 
in a Luttinger liquids [§^,0 and some related conjectures on non-Fermi liquid aspects of 
high Tc superconductivity, we elaborate on the potential impact of the results and discuss 
their range of validity. An appendix contains many of the technical details. 

II. HIGHLIGHTS OF A FERMI LIQUID 

The most succesful description of interacting degenerate Fermi systems in the normal 
state, i.e, non-superconducting or superfluid is Fermi liquid theory. Landau's original for- 
mulation, largely phenomenological, has as basic hypothesis a one-to-one correspondence 
between the eigenstates of the interacting and non-interacting system. In this formulation 
the quasipartides are obtained from the single particle states of the non-interacting system 
via an adiabatic switching-on of the perturbation PU|-P[], hence for this picture to remain 
valid the interactions should not lead to phase transitions. The Landau quasiparticle concept 
is appropriate for excitations very near the Fermi surface, since for short range interactions 



the lifetime of these quasipartides is r ~ |/c — /ci?| p9|-p|, which after including screening 



effects also holds in the case of Coulomb interactions [43|. For quasipartides near the Fermi 



surface the adiabatic hypothesis is reasonable [^,|^ and Landau's phenomenological theory 



is applicable to study transport phenomena of low energy excitations ^T| . A more modern 
and consistent description of Fermi liquid theory is based on renormalization group ideas | H 
which describe the low energy physics near the Fermi surface as fixed points of the renor- 
malization group, and Fermi-liquid interactions are those associated with marginal operators 
near this fixed point. This formulation of Fermi liquid theory reveals that the physics near 
the Fermi surface is very similar to that of critical phenomena and is completely determined 
by the gapless excitations associated with the formation of particle-hole states near the Fermi 
surface @ . 

Fermi liquid theory is the starting point of a consistent study of the properties of degen- 
erate, interacting Fermi systems, in particular BCS superconductivity can be understood as 



a result of the instability of the Fermi liquid towards pairing attractive interactions [^ 



It has been adapted to study dense nuclear matter [^] and recently the renormalization 
group description of Fermi liquids has been extended to the interactions in field theory at 
nuclear matter density [0. The conclusion of these studies is that Fermi liquid theory is 
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the effective low energy theory for excitations near the Fermi surface. The field theoretical 
approach to Fermi liquids quantifies the main concepts of Fermi liquid theory in terms of 
the spectral density for the quasiparticles near the Fermi surface: a Breit-Wigner form with 
a small width (for k ~ kp) and a finite residue at the "quasiparticle pole", Zkp. The 
quasiparticle distribution function is distorted from the original Fermi-Dirac step function 
(at zero temperature) but a "jump discontinuity" remains at the Fermi surface which is 
determined by the wave function renormalization constant Z^^, (see ^^ for an explicit 



calculation in the electron gas). In a well defined sense, this jump discontinuity is associated 
with an "order parameter" for a Fermi liquid [40,^|: a normal Fermi liquid corresponds to 



Zkp 7^ while for a non-Fermi liquid Z^^ = 0. Transport properties of a degenerate Fermi 
gas depend on the renormalization of the (quasi) particles, for example the coefficient of the 
linear power of temperature in the electronic specific heat is proportional to Z^p- 

There are some notable examples of the breakdown of Fermi liquid theory such as the 
Kondo model of electrons interacting with magnetic impurities, and "quasi" one dimensional 
metals which provide a novel type of behavior for correlated electron systems: the Luttinger 



liquid (see |3^,5D] and references therein). This novel behavior of a correlated degenerate 
electron system is characterized by a vanishing jump at the Fermi surface and power law 
correlations akin to those found in critical phenomena ||39(| . 

Another system in which the Fermi liquid description breaks down is that of non- 
relativistic electrons interacting via the exchange of "magnetic" (transverse) gauge bosons 

or 



50|. Recent conjectures suggest that this type of non- Fermi liquid behavior [JH 
"Luttinger-liquid" behavior could explain the unusual properties of the normal phase of 
high temperature superconductors P9| . 

Our study of the normal state of cold-dense QCD (and of QED) reveals a striking re- 
semblance with Luttinger liquid behavior, in particular the renormalization group improved 
quark propagator that we find is remarkably similar to that proposed in |58.|^ to describe 
the normal phase of high temperature superconductors. 

III. PRELIMINARIES: QUASI PARTICLE DISTRIBUTION FUNCTION AND 

EQUATIONS OF MOTION 

In this section we obtain the general expressions for the quasiparticle and quasi anti- 
particle distribution functions to make contact with Fermi liquid theory. The main aspect 
of the distribution function is that the existence of a jump discontinuity of the quasiparticle 
distribution at the Fermi momentum is the signal of Fermi liquid behavior. Furthermore 
we obtain the effective equation of motion for the Dirac field in the medium to extract 
non-equilibrium aspects. Typically relaxation es studied by extracting the damping rate 
which is obtained from the imaginary part of the self-energy evaluated on the mass shell 
of the Fermion and describes exponential relaxation. However, from the results of p^^7| 



it is found that the damping rate vanishes for quasiparticles with the Fermi momentum. 
However the form of the Fermion self-energy (see next section) strongly suggests the build- 
up of logarithmic infrared singularities just like in critical phenomena with the potential 
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for suinming up to a power law relaxation with anomalous dimensions. The solution of 
the real-time equation of motion resummed via the dynamical renormalization group will 
confirm power law relaxation with anomalous dimension for single quasiparticles with the 
Fermi momentum. Thus in order to extract the correct relaxation behavior we must study 
the real-time evolution of the amplitude of the quasiparticle state which is obtained from 
the equation of motion. 

A. Quasi particle distribution function 



The spatial Fourier transform of the Dirac field operator at any given time t can be 
written in the form 



mt) = y: [fe^wf/^^n^) + diit)v^^^\-k) 



(3.1: 



with U^'^'{k);V^'^'{k) the usual free particle Dirac spinors. We refer to the time depen- 
dent operators b^sit);dl (t) as the annihilation and creation of quasi-particles and quasi- 
antiparticles respectively. Within the spirit of Fermi liquid theory, upon adiabatically 
switching-on of the interaction these operators interpolate between the free (bare) parti- 
cles and the dressed (quasi) particles and antiparticles respectively. We define the average 
number of quasi-particles and quasi-antiparticles as 



rii 



ni 



^ s 



(3.2) 
(3.3) 



Where the expectation value is in the exact ground state, which is obtained from the 
unperturbed ground state by adiabatically switching on the perturbation from time t — *> — 00 
to t = 0. Using the properties of the usual spinor wave functions [/, V and the results of the 
appendix, it is straightforward to find 
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4cUr 
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^ |' 7o(^ + m)7o 
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(3.4) 
(3.5) 



We are considering the situation of fi ^ m so that we will neglect the current quark 
masses and consider the quarks as massless. Using the spectral representation of the propa- 
gators for Dirac fields given in the appendix, we find for massless fermions 



ni 



rii 



- I dqoTr \V+{k)pf{qo, k)\ Nf{qo) 



- I dqoTr V-{k)pf{qo,k) 



Nfiqo)) 



(3.6) 
(3.7) 



With V±{k) given in the appendix. Thus to obtain the quasi (anti) particle distribution 
functions we need to obtain the fermion spectral density p(go, k). 

The emission and absorption of hard gluons (photons) with momenta q > fi will affect 
the distribution functions for fermions from deep within the fermi sea up to the excitations 
near the fermi surface. However, because of Pauli blocking of states below the Fermi surface, 
the emission and absorption of soft gluons (photons) with q << fi will only affect the 
distribution function of particles near the Fermi surface but not those deep within the Fermi 
sea. Soft gluons (photons) are sensitive to screening corrections and their propagator must 
include screening arising from quark loops. The leading order correction is given by the 
resummation of hard dense loops (HDL) corresponding to quark intermediate states with 
momenta near the Fermi surface ||3^- p7|^ . The spectral density in both cases, with hard 
and soft gluon exchange, will be studied in detail in the next section. 

B. The equation of motion 

We can treat equilibrium and non-equilibrium aspects of cold and dense QED and QCD 
by studying the real time equation of motion of a fermion condensate induced by an external 
source. The equilibrium aspects studied here can all be addressed by obtaining the fermion 
propagator and the fermionic spectral density, while the equation of motion allows us to 
study the real time relaxation of an initially prepared condensate as an initial value problem 
prepared by using a suitable source term. The advantage of studying the equation of motion 
in real time is that it will reveal the relaxation of fermions directly in real time. Since 
the leading order corrections in the HDL approximation are similar for QED and QCD 



t ^ and the gluon (photon) polarization is completely determined by one fermion 
loop, we describe the necessary steps in QED, the final form for the spectral densities and 
relevant quantities for QCD can be obtained by simply accounting for the proper color and 
flavor factors. 

The QED Lagrangian density in Coulomb gauge is given by (see [^] for a similar context) 



C = ^ (i^~ g-f^A^ + g^ ■ At) ^ + ^77 + 77^ 



+ 2 



{d^ATf + iVA')' 



(S.t 



where the Grassmann valued source terms were introduced to obtain the effective Dirac 
equation in the medium by analyzing the linear response to these sources. 

Following the steps detailed in |^, we find the Dirac equation for the spatial Fourier 
transform of the expectation value in the massless case to be given by 

M7o^ -1-kj ^{k, t) + y_^ dt' S(A;, t - t') ^{k, t') = -r]{k, t), (3.9) 

where S(A;, t — t') is the retarded fermion self-energy given by the sum of the transverse and 
longitudinal contributions 
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S(A?, t-t') = St(A?, t-t') + Sl(A?, t - t') 
Using the results of the appendix we find 



T.T{k,t-t')=ig^ 



(27r) 



V^Av){l'[{^S^^it,m-^g^t{t,t')) 



(3.10) 

{^S^{t,t')){-^G<^^{t,t'))\^^] 
(3.11) 



2 / C?3? 



(3.12) 

with p = k — q. The propagators are written in terms of their spectral representation and 
using the results in the appendix it becomes clear that the retarded self-energy has the 
following causal structure 



S(A?, t-t') = T:(k, t - t')Q{t - f) 



(3.13) 



Introducing the Fourier transforms in time for the expectation value of the Dirac field, 
the source and the self-energy, and the Fourier representation of 6(t — t') we find that the 
equation of motion in terms of the space-time Fourier transforms becomes 



7°cj — 7 ■ /c + S(cJ, k) ip{uj, k) = —rj(uj, k) 
Introducing the dispersive representation for the retarded self-energy 

S^, k) = - dqo -— 

71 J qo — LO — iU+ 

a straightforward calculation reveals that 

E'-(A?, t-t') = - f dgoe-''^°(*-*')lmE(go, k) 
For massless fermions we write 

S(a;, k) = 7*'So(c<J, A;) — 7 ■ A;Si(cj, k) 
Hence the solution of the equation of motion is given by 

i^iuJ^k) = SR{uj,k)fi{u,k) 
with 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



10 



Sr{u, k) 



V4k)S4uj, k) + V+{k)S+{uj, k) 



S^ioj.k) = uj — k + T,_{uj,k) 
S+{uj,k) = uj + k + T,+ {uj,k) 

S_(ti;, k) = So(^5 k) — Si(ti;, k) 
I]_|_(u;, k) = T,q{u, k) + T,i{u, k) 



-1 
-1 



(3.19) 
(3.20) 
(3.21) 

(3.22) 
(3.23) 



with V±{k) given in the appendix. The fermion spectral density is obtained from the imag- 
inary part of the retarded propagator and is given by 



Pf{u, k) = - \P.{k)p.{uj, k) + P+(A;)p+(^, k) 



p-{uj,k) 



ImS_(a;, k) 



TT 



p+{uj,k) = - 



io — k + ReT,_{uj,k) + ImS_(ci;, A;) 
ImS+(ti;, k) 



vr 



cij + A; + ReS+(ci;, A;) + ImS+(c(j, A;) 



(3.24) 
(3.25) 



We consider two cases separately to obtain the spectral densities: i) the gluon (photon) 
line in the fermion self energy carries hard spatial momentum p > p. In this case one 
(bare) gluon exchange gives the leading order correction to the quark self-energy in the HDL 
approximation pij-jT^ for large chemical potential, ii) the gluon (photon) line carries soft 
spatial momentum p <^ p in which case the gluon (photon) propagator must be dressed by 
HDL fermion loops |34-^. The contribution from hard gluon exchange to the self-energy of 
low momentum fermions deep within the Fermi sea, i.e, with k « gp, must be treated non- 
perturbatively 



l2|,p4H37[| resulting in a modified dispersion relation and their description as 
quasiparticles. For fermions near the Fermi surface (and weak coupling) the HDL corrections 
from hard gluon exchange are perturbative. Thus hard gauge fields will modify the fermion 
propagators for all fermion states in the Fermi sea. On the other hand the emission (and 
absorption) of soft gluons (photons) with q < gp will only affect fermionic states near the 
Fermi surface. We now study each case in detail. 



IV. EQUILIBRIUM ASPECTS 
A. Hard gauge fields (HDL): Fermi-liquid behavior 

We are interested on particles near the Fermi surface, therefore we will focus on the 
spectral density p-{qo, k) and we will neglect the anti-particles, for which there is no Fermi 
surface. 

For hard gauge fields with loop momentum q > p we can use the free field propagators 
for transverse and "longitudinal" gauge fields given in the appendix. We find that the 
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instantaneous Coulomb interaction leads to a local contribution to the self-energy which is 
subleading for small fermion momentum [^,^]. In the HDL limit we find 



p-iqo,k) = Z+6{qo - uj+{k)) + Z_5{qo + uj^{k)) + pc{qo,k) 



(4.1) 



with uj+{k) , uj-{k) the fermion and plasmino quasiparticle poles with residues Z± respec- 
tively (see [Q for their complete expressions). The continuum spectral density has support 
below the light cone and is given by 



Pciqo, k) 



ImS_(ci;, k) 



vr 



~ — * Zi — — * 

u; — fc + ReS_(co', fc) + ImS_(ti;,A;) 



lmI]_(a;,^) = 7r^(l-^)e(A;2 



ReS_(u;, k) 



Mj 



and satisfies the sum rule 



_Ml 
2k 

87r2 



iU 



In 



uj + k 



LU 



+ 2 



Z_ + 



■iqo, k) 



(4.2) 

(4.3) 
(4.4) 
(4.5) 

(4.6) 



ni 



riT 



The spectral densities p±{qo, k) are related by |Q 

P+{qo,k) = p^{-qo,k) (4.7) 

and upon using this relation, the results of the previous section and the appendix we find 

Xqo,k)N{qo) (4.8) 

loP-{Qo,k)N{qo) (4.9) 

N{qo) = e{p-qo) ; N{qo) = &{-p - qo) (4.10) 

and the fermion number density is given by 

d^k 



AT 



(2vr)= 



dqop-iqo,k) 



(4.11) 



this expression relates the chemical potential to the fermion number density. We now have 
the tools to understand the change in the Fermi sea in the HDL limit. 

Since N{qo) = Q{p — q^), and uj^{k) > k |3^, from the expression for the particle 
distribution function given by ( [4.8|) , we see that for values of k such that L(j+(fc) < p the 
region of support of the spectral density P-{qo, k) is contained in the interval — oo < qo < p, 
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therefore for these values of k the distribution function is n^ = 1. The value of the momentum 
at which the frequency of the fermion quasiparticle uj+{k) = /i is the limiting value for 
which the quasiparticle pole is in the interval of support of N{qo), thus defining the Fermi 
momentum kp by u^lkp) = /i. 

We can obtain an estimate for the value of the Fermi momentum in the HDL limit and 
weak coupling for which fi ^ Mf using the large k limit of the quasiparticle pole 



uj+{k) ^ k + —f- k:^ Mf 



from which we obtain 



kp K. ^ 



1- 



^vr" 



(4.12) 



(4.13) 



For k > kp the quasiparticle pole is outside of the region of support of iV(go) but the 
spectral density still has a contribution inside this region, given by the plasmino pole at 
—uj_{k) and the Landau damping continuum. Hence nfc>fc^ 7^ as befits an interacting 
Fermi system. Thus whereas for k < kp the distribution function ra^ = 1, for k > kp it is 



< ng < 1. 



This analysis leads to the following result 



n^ = Z_{k) + Z+{k)Q{kp -k) + 



Pc{qo, k) 



1 k = kp — e 
l-Z+{kp) k 



kp + e 



(4.14) 



We can estimate the discontinuity or "jump" at the Fermi surface in the HDL limit by 
using the sum rule ([4.6|) above and the large k limit of the quasiparticle residue ||33| 



1 - Z+{k ^ kp) 



At the Fermi momentum we find 



An) 



2 In 



47r 



-1 



nk=kF = N{qo = fi)Z+{kp) + f ^j 




(4.15) 



(4.16) 



where we have used N{qQ = n) = 1/2. Thus to this order, the Fermi surface is sharp with a 
jump discontinuity at k = kp given by nk=kp-e — nk=kp+e = Z+{kp). The sharpness of the 
Fermi surface in the HDL approximation is a consequence that to this order the quasiparticle 
excitations have an infinite lifetime. For momenta k ^ fi the HDL approximation is not truly 
justified since in this limit k ^ g^x and the contribution from the hard momentum region 
of the two particle cut is of 0{g'^fi) which becomes comparable to that of Landau damping 
of 0{g'^fi^/k). This fact notwithstanding, the main purpose of this analysis is to reveal that 
the contribution from the hard loop momentum region to lowest order in the resummation 
program leads to a description consistent with Fermi liquid theory with a discontinuity or 
"jump" at the Fermi surface determined by the residue of the quasiparticle pole which to this 
order is non-vanishing. As we will study in detail below, the contribution from soft gluon 
exchange, which requires screening corrections invalidates the Fermi liquid description. 



13 



B. Soft gauge fields: Non-Fermi Liquid 



We now focus on the damping effects on quasiparticles near the Fermi surface, i.e. k ~ 

kp ~ /i. Damping of these excitations results from the emission and absorption of soft gluons 

which require the gluon (photon) propagators in the fermion self-energy to be HDL resummed 

3^-0 . These dressed gauge propagators can be handily included in the calculation of the 



self-energy by writing the Wightman functions for the gauge fields in terms of their spectral 
representation described in the appendix. Since the fermion momentum is A; ~ fci? ~ /i and 
the exchanged gluon is soft with g -C /i, the fermion line in the self-energy does not require 
HDL resummation and can be taken to be a bare fermion propagator. 

A straightforward calculation using the free fermion propagators for the internal fermion 
line and the HDL resummed longitudinal and transverse gauge field propagators in terms of 
their spectral densities leads to the fermion self energy in a dispersive representation as in 
i ^JBi ) with 

ImE(go, ^) = ^ / 7^ / ^Po{[e(po) - e(/i - q)] 6{qo ~ Po - q) [pt{po, p)Vi^ {p)YV.{^W + 



PL(po,p)7°^-(q)7°] + Q{Po)5{qo+Po + q) 

pl{p,.p)v'4{y>WvM)i'+Pl{p,.p)i'vM)i']} (4.17) 



where p = k — q and we have neglected the instantaneous Coulomb interaction and pt,l are 
given in the appendix. 

The dispersive representation ( p.l5| ) makes clear that the largest contribution to the 
self-energy for cj ^ /i is determined by the behavior of ImS(go, k) for q^ ^ u ^ fi. 

The contribution from soft gluons to the self energy of fermions near the Fermi surface 
k ^ fi can be extracted easily by first relabelling the integration momenta p ^^ k — qso that 
g is the momentum of the (soft) gluon line with g -C /i and \k — q\ ^ k — qcos{6) with 6 the 
angle between k and q. In this limit we find 



ImS, 



d^q 



.(go, k) = ng^ j — -^ j dpo { [pt{po, g)(l - cos^(^)) + pl{po, q) 

(e(po) - e(/i - A; + g cos(^))) <5(go - po - k + q cos(^)) + 
pApo, g)(l + cos\9))e{po)S{qo + Po + k- gcos(^))} 



(4.18) 



ImE+(go,A;) = ng^ 



d^q 



c^Po {pt (j5o,g)(l + cos^(6')); 



(e(po) - e(p - /c + qcos{9))) 6{qo - po - k + q cos(e)) + 
Pt{po, g)(l - cos2(0)) + pl{po, g)] Q{po)S{qo + po + k - gcos(^))} (4.19) 

For particles near the Fermi surface the quark propagator has poles near u ^ k ^ p, while 
for antiparticles there is no Fermi surface and the poles are at cj ~ —k. Therefore for particles 
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the self energy is determined by IniS_|_(go, fc) for go ~ ^ through the dispersion relation. 
Hence for particles near the Fermi surface the important region is go ~ ^ ~ ^ which implies 
that the argument of 6{qo — po — k + qcos{9)) has support for po ^ qcos{9) corresponding 
to the region of Landau damping of the gluon (photon) propagator (see appendix). On the 
other hand ioi qo ^ k ^ u ^ fj, the 6{qo + po + k — qcos{9)) has support in the region of 
Po ~ —2fi < and this contribution is cancelled by the 0(po). Thus the contribution to the 
self-energy of particles near the Fermi momentum from soft gluon exchange is completely 
determined by the first delta function 



ImS_(go,A;) ^ ng^ 



d^q 

(2^ 



{ [pt(po, g)(i - cos^(^)) + pl(po, q) 



[©(go -k + qcos{9)) -Q{fi-k + gcos(^))]} 



po=qo-k+qcosl 



(4.20) 



which obviously vanishes at go = //. Thus we see that the imaginary part of the self energy 
i.e. the damping rate for fermionic excitations vanishes at the Fermi surface |^B|,^. For 
antiparticles the propagator has poles for uj ^ —k, for hard momentum k ^ fi and go ~ —k 
the 6{qo — po — k + qcos{9)) has support for po ~ —'^k ^ — 2/i, i.e, in the hard region, while 
5(go + Po + k — qcos{9)) has support for po = —go — k + qcos{9) ~ gcos(6') which is in the 
Landau damping region of the spectral density. Hence 



ImS+(go, A;) ^ irg'^ 



d^q 

^2^ 



{ fpr(Po, g)(i - cos^(^)) + Pl(Po, q)] 0(Po)} 



-qQ—k+qcos(9) 

(4.21) 



The case of antiparticles has been studied in P7| , |35| with the result that the self-energy 
is analytic near the Fermi surface and does not lead to novel phenomena. Hence in what 



follows we will ignore the case of antiparticles and refer the reader to |37,35| for more details 
on this case. 

Since the imaginary part of the self-energy for particles vanishes at the Fermi surface, we 
expand near go ~ yU as follows 



6(go — k + qcos{9)) — 6(/U — k + qcos{9)) = qoS{qcos{9) — (A; — /i)) + 
with 

•fo = (go - P) 
and we are led to the following expression 



ImS_(go,/i;) ^ -— g~o / 

47r J\k-^l\ 

For go —^ we can approximate 



PT(go, g) 1 5 — + PL(go, g) 



dq 



(4.22) 



(4.23) 



(4.24) 
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PT{qo,q) 
pL{qo,q) 



M' 



2/90 



)eif - qV) 



g2 + 4M2(|^ 



+ 



1 



2M2(2o) 



[g2 + 4M2]^ 

47]-2 



(4.25) 

(4.26) 

(4.27) 



The region in which dynamical screening via Landau damping is effective is determined 
by 5', Iffol < M, hence the vahdity of the approximation invoked above rehes on |(fo| < M. 

We note that whereas pt has a strong infrared singularity for g ^ when go ^ 0, 
Debye screening cuts off the infrared behavior of pi which leads to the conclusion that 
the contribution from the longitudinal gluons (photons) is oc go^ [131 • -^^^ excitations at the 
Fermi surface with k = p the contribution to the imaginary part from the transverse photons 
can be computed straighforwardly and it yields 



ImSL")(go,fc) 



24:71 



ml 



(4.28) 



The contribution from the longitudinal gluons (photons) is clearly of 0{g^qQ^) and of the 
same order as the contribution from the hard gluon (photon) region [R^. This is similar to 



the case of non-relativistic electrons interacting via the screened Coulomb interaction [^ 
The term proportional to k = k — p in ( [4.24| ) can be easily shown to yield a contribution of 
Oig^{kyM)\C[;/M\^/^). 

Therefore we conclude that 



Imi;_(go,/c) 



247r 



^ 2 



\qo 



2^0 



^^^'it 




UJ 



M 



(4.29) 



with calculable coefficients A, B. Our focus is to understand the quasiparticle excitations 
and their distribution function very near the Fermi surface, in particular the discontinuity 
of the distribution function at kp- In order to do this we need only consider k ^ M with 
|ci)| -C M therefore we will keep only the first (leading) contribution in ( |4.29|) . 

We can now obtain the self-energy via the dispersion relation ( 3.15|) ioi uo k. p by using 
the first term of ([4.29|) and integrating within a region of width ^ M around p since this is 
the region in which Landau damping is effective for dynamical screening and the region that 
yields the leading infrared contribution. 

We finally obtain for fc, go ~ /i 



S_(ci;, k = p) 



uo 



»' ak 


1 


M 


+ ln 


"cD + zO+" 
M 


247r2^ y 


g' ~, 


UJ 

M 


+ ^24^^" 


uo — p 













(4.30) 
(4.31) 
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We finally obtain the retarded propagator and spectral densities for particles near the 
Fermi surface by using the relations (|3.19| - |3.25|) and extracting the term proportional to 
V-{k). 

It is convenient to introduce k = k — fi^k — kp (since from the analysis of the previous 
section kp = fj,{l — 0{g'^))), for u,k ^ fi the inverse propagators for particles can be written 
as 



S_'^{uj, k) ^ LU — k + T,^{lj, k = fi) 



(4.32) 



the form of the real part of the self energy oc cDln ItD/Ml strongly suggests a wave function 
renormalization. Such an interpretation is obscured by the term k in the propagator, how- 
ever, at the Fermi surface with A; = 0, the logarithmic singularities of infrared origin are 
reminiscent of those of a critical theory and the propagator is similar to that obtained in 
perturbation theory in a critical theory at the upper critical dimensionality. This interpre- 
tation is validated by the current ideas on Fermi liquid theory based on the renormalization 
group, which describes the effective theory near the Fermi surface as a critical theory with 
marginal Fermi liquid interactions [0 . 

Since at the Fermi surface {k = 0) the inverse propagator is proportional to uj the wave 
function renormalization constant at the (quasi) particle pole for excitations near the Fermi 
surface would be given by 



Z(/i) 



1 + — ReE_(t^, k)\^=k=fi 



1 -1 



•"^ 




1 




1- 


- 7^2 In ^ " 



(4.33) 



This wave function renormalization or residue at the (quasi) particle pole for excitations 
near the Fermi surface is precisely the quantity that determines the "jump" of the (quasi) 
particle distribution function at the Fermi surface as determined by ( [4.16|) . However, the 
logarithmic singularities manifest in ( [4.30|) would lead to the conclusion that Z{fi) = and 
that the Fermi surface "vanishes" . Such conclusion has also been obtained in non-relativistic 
systems with magnetic interactions p8H30[]. 

A similar situation arises in a critical theory, for example in the Euclidean formulation 
of a critical scalar theory with quartic interaction in four dimensions the inverse propagator 
for small Euclidean four momentum K is given by 



G-\K)^K' 1-X'c\yi{K'/k' 



(4.34) 



with A the quartic coupling, c a combinatoric constant and k a renormalization scale. Again 
the wave function renormalization or residue at the "pole" K^ = 0, 



d 



dK- 



■G-\K)\k2=o 



-1 



(4.35) 



vanishes. The vanishing of the wave function renormalization on the particle mass shell in a 
critical theory indicates that the propagator acquires an anomalous scaling dimension. The 
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logarithmic singularities are resummed via the renormalization group leading to the following 
improved propagator 



G-^UK) = K' 



K' 



K^ 



V 



\\ 



(4.36) 



Therefore we interpret the vanishing of the wave function renormalization for quasipar- 
ticles at the Fermi surface as an indication of the build-up of an anomalous dimension as a 
consequence of the infrared singularities, which are only dynamically screened via Landau 
damping. To make this interpretation explicit we now proceed to obtain a renormalization 
group improved quasiparticle propagator, focusing solely on the particle excitations near the 
Fermi surface. 



C. Euclidean Renormalization Group 



i. fc = o 



In order to obtain a renormalization group resummation of the infrared singularities we 
must first perform an analitical continuation to Euclidean space. This is accomplished by 
the following analytical continuation 

cD + iO+ = iK (4.37) 

with K taken to be a real variable. From ( p.20|) the particle propagator now reads 



S-{u:,k)\^r^o^k=o 



K 



T{K) 



T{K) = K^ 



A In 



M2 



A 



247r2 



(4.38) 
(4.39) 



Obviously T{K) has the same form as the inverse propagator of a scalar theory with 
infrared corrections typical of a critical theory as discussed above. The physics near the 
Fermi surface requires K <^ M, which obviously leads to a breakdown of the perturbative 
expansion. Just as in a critical theory taking K —^ keeping M fixed is the same as keeping 
K fixed and taking M -^ oo, i.e, interpreting M as an ultraviolet cutoff and taking the limit 
of large cutoff at a fixed transfered momentum. 

The renormalization group improvement proceeds in the same manner as in the scalar 
theory, we first introduce a wavefunction renormalization constant that will absorb the cutoff 
dependence at a given renormalization scale k, and a renormalized vertex function 



Tr{K,k) = Z[k,M]T{K,M) 
Z[k] = 1 + Aln ^^ + 



M2 



(4.40) 
(4.41) 



The independence of the bare vertex upon the renormahzation scale leads to the renor- 
malization group equation 



d 



r^(ir,«:) = o 



KdZ 

Using the fact that the vertex has scaling dimension 2 i.e, 



Tr(K.k) 



"'" ' t) 



(4.42) 
(4.43) 

(4.44) 



with $ a dimensionless function of its argument, the renormahzation group equation for the 
scaling function $ can be solved straightforwardly leading to 



rR{K,K) = ^ 



K 



2~ri 



$fll 



(4.45) 



We can now obtain a renormahzation group improved vertex function as follows. Re- 
quiring that at the scale K = M (i.e, at the scale of the cutoff where perturbation theory is 
valid) the "bare" vertex obeys 



T{K = M,M) = M' = Z-'[k,M]Vr[K = M,K 



(4.46) 



fixes the renormalization group resummed wave function renormahzation in terms of the solu- 
tion of the renormalization group equation ( [4.45| ) at the scale K = M. The renormalization 
group improved "bare vertex" T{K, M) = Z^^[k, M]rfj[K, k] is valid for K <^ M and leads 
to the following resummed particle propagator valid ioT k = kp, uj < M and weak coupling 
A« 1, 



^-(^,^)l 



u^O,k=0 



I 

'm 



M2 



(4.47) 



The resummed spectral density for particles p-{u!,k) = lm.S-{uj,k)/7c is obtained by 
performing the analitic continuation K -^ —iuj + O"*" leading to the following expression near 
the Fermi surface 



p-{u,k)\ 



LUKifl,k=fl 



sin[7rA] 



TTlUJl 



LU 



M 



2X 



u 



u 



f^ 



A 



9 



247r2 



(4.48) 



It is straightforward to check that this spectral density becomes a delta function in the 
limit A ^ by integrating it within a small region around the Fermi surface. This spectral 
density is remarkably similar to those found in non-fermi liquid systems such as Luttinger 
liquids and have been experimentally measured in condensed matter systems via the X-ray 
edge singularities at the Fermi surface of some metals |H2 . 
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Using the renormalization group improved spectral density ( [4 .481 ) we can now obtain the 
value of the particle distribution function at the Fermi surface by restricting the integral in 
( [4.8| ) to a region of width M near the Fermi surface 



n 



k=kf 



" du Qi-uf''^''^^ 



-M 



TXU 



U 



M 



2A 



N, 



sin[7rA] 



k=kf 



71 



A 



(4.49) 



with A^^ the Fermi-Dirac distribution function and A^^^^ 
( [4.49| ) obviously attains the free field limit for A ^ 0. 



1/2. The distribution function 



2. k^O 



The renormalization group method presented above for the case oi k = can be straight- 
forwardly extended to k ^ hj recognizing that the propagator for A; 7^ is 



T{K) = K^ 



A In 



M2 



+ iKk. 



(4.50) 



similar to that obtained in a scalar theory in the large transfered momentum limit but with a 
"mass" term. This similarity suggests that the term ikK can be treated just as a mass term 
in the renormalization program of a scalar theory. Therefore, along with the wave function 
renormalization ( |4.40| ) we also introduce the following multiplicative renormalization 



kR = Z[K,M]k 



(4.51) 



since the are no infrared divergences associated with k then Z[k,, M] is chosen in perturbation 
theory to be that of the k = case discussed above and is independent of k. 
The renormalization group equation now becomes 



d ~ d 

OK dkn 



TniK, kn, k) = 



V 



7 



KdZ _ 
Z dn 
K, dkn 
kn Ok 



2A 



2A 



using the fact that the vertex function has dimension 2 we write 



Tr{K, kn, k) 



\ K K 



Finally introducing the following variables and ansatze 



K _ 

K 

$(e-*, 



kn 



K 



e*(''-2)i7[t, 



(4.52) 
(4.53) 
(4.54) 

(4.55) 



(4.56) 
(4.57) 
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the function H\t, r] obeys the simple renormahzation group equation 



d 



d 



dt dr 



H[t,r] = 
7 = 7-1 = 2A-1 
The solution of this equation is simple 

H[t,r] = H[f{t)] 
with f (t) the solution of the ordinary differential equation 

dv 

= -7f =^ fit) = f(0)e"^* 
dt 



(4.58) 
(4.59) 

(4.60) 

(4.61) 



Therefore we find that the renormahzation group improved vertex function is given by 



TR{K,~k,K) = K' 



'K^V^ 



H 






(4.62) 



where we have explicitly used the multiplicative renormahzation of k and introduced the 
integration constant Z that only depends on n and M to be determined later. 

Comparison with the solution ([4.45|) for A; = reveals that if [0] = $(1). As in the case 
fc = we request that the "bare" vertex coincides with the free field expression at the cutoff 
scale K = M, i.e, 



T[M,k,M] =M^ 



l + t 



M 



Z-'[K,M]Tn[M,k,K\ 



(4.63) 



this condition and the identification H[0] = $(1) determine both Z[k,,M] and Z[k,,M] by 
comparing the powers of k . Since if is a function of the scaling variable ip = kZ{K / kY^ the 
functional form of H[(p\ is uniquely determined for k <^ M. We are thus led to the following 
unique form of the renormahzation group resummed vertex function for fc <^ M ; A ^ 1 



T{K,k) = K^ 



M2 



-A 



1 + 



ikK 



K^m 



-A 



(4.64) 



This discussion reveals another manifestation of the critical nature of the theory near the 
Fermi surface: the renormahzation group improved vertex function is of the scaling form 



T(K,k) = T{K,0)D 



k 



A = 1-2A 



(4.65) 



Using the relation between the vertex function and the Euclidean particle propagator 
given by ( |4.38| ) we find from ( |4.64| ) 
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S-{iJ^ k) 



uj^fi,k'^fi 



'iK m 



-2A 



K=-iCo + 0+. 



k 



(4.66) 



we can now obtain the spectral density near the Fermi surface by following the steps 
detailed in the case fc = 0, leading to 



P-{uJ,k)\^^^^k 



'll,KP^ll 



sin[7rA] _ 



vr 



UJ 



M 



-2A 



/ 


M 


-2A 


— A;cos[7rA] 


2 ,~ ^ s2 

+ f A;sin[7rA]j 



(4.67) 



which vanishes identically at cj = 0. This spectral density is displayed as a function of u for 
several values of k and A in figures (1) and (2). 



3. Quasiparticles 

Figures (1) and (2) clearly reveal that in weak coupling the spectral density features a 
narrow resonance for fc 7^ 0, the position of which is obtained by the vanishing of the first 
term in the denominator in ( [4.67|) which determines the quasipartide dispersion relation. It 
is given by 



UpCk) = sign(Ai) [|A;|M-2^cos[7rA]] '"'" 
We note that the group velocity of the quasiparticles near the Fermi surface 



Vg{k) = M-'-" cos[7rA] 



~ 2A 



(4.68) 



(4.69) 



vanishes as k ^ kp. We interpret this novel phenomenon in terms of a collective backflow 
that surrounds the quasipartide. 

Near the position of the resonance, i.e, for Co ^ Up the spectral density can be approxi- 
mated by a Breit-Wigner form 



p_{uj, k)\cjr^cjp,k^f, = Zp[k\- 



m 



n 



{u-Up{k))'' + T^{k) 



Zp[k] 



M 



2\ 



(1-2A) 



r(A;) = Zp[A;]|A;|sin[7rA] 



(4.70) 

(4.71) 
(4.72) 



Therefore the residue of the "quasipartide pole" and the "quasipartide width" vanish 
near the fermi surface as 
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Zp[k] oc |A; — kpl i- 
r(^) oc \k-kF\~ 



2A_ 

2A 



(4.73) 
(4.74) 



It is straightforward to confirm that the quark propagator ( 4.66| ) has a complex pole 
with the real and imaginary parts given by ( |4.68|) and ( [4.72|) respectively in the narrow 
width approximation T{k)/up{k) <C 1. For this purpose write for the complex pole 



Kp = —iujp — r 



(4.75) 



with r > corresponding to damping (this is confirmed a posteriori from the solution). In 
the narrow width approximation we can replace 



{ijjp — iV 



,1-2A 



r 



{cjpy-'^\i-{i-2\)i^ + 



UJr. 



(4.76) 



requiring the vanishing of the real and imaginary parts of the denominator of ( |4.66| ) we find 
the position of the quasiparticle pole and its width given by ( [4.68| ) and (|4.72|) respectively. 



4-- Requiem to the Fermi Liquid: the jum,p of the distribution function at kp vanishes 

We can now study the behavior of the distribution function near the Fermi surface by 
using the spectral density (|4.67| ) in the expression (|47^ ). Since the domain of validity of the 
approximations invoked to obtain the spectral density is such that \uj\ , \k\ < M the integral 
in the variable uj must be cutoff at the scale M. Hence we find 



ni 



Uj 



sin[7rA] 



27rA 



-M 



An{k) 

A detailed calculation reveals that 
An{k) ■■ 



p^{uj,k) - p^{uj,0) 



(4.77) 
(4.78) 
(4.79) 



k ^ {-k/MY sin [{n + 2)7rA] 



ttM 

~k 
"ttM 



n=0 



n + l-2{n + 2)X 



sin27rA k sIuSttA k"^ sin47rA ^o/rsx 
1-4A ~ M(2 - 6A) ^ M2(3-8A) ^ ^ ^ 



(4.80) 



leading to the following form of the single quasiparticle distribution near the Fermi momen- 
tum 



n 



koikf 



sin[7rA] 
27rA 



k 

7rM(l-4A) 



+ 0(P 



(4.81) 
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This is one of the important results of this work: the exchange of magnetic gluons, which 
are only dynamically screened by Landau damping are responsible for the vanishing of the 
discontinuity of the distribution function at the Fermi surface. This is clearly a consequence 
of the result that the wave function renormalization or residue at the quasiparticle pole 
vanishes as a power law at the Fermi surface. In Fermi liquid theory, the discontinuity at the 
Fermi surface enters in all thermodynamic response functions, for example in the coefficient 
of the linear power of temperature in the specific heat ||4l| , |42| , thus the vanishing of the jump 



discontinuity will likely result in an anomalous specific heat. 

V. NON-EQUILIBRIUM ASPECTS: SINGLE QUASIPARTICLE RELAXATION 

A. Dynamical Renormalization Group 

The equation of motion ( p.9| ) can be written in the form 

with 



no^ -1-k] ^{k, t)+ I __ dt' Z'-{k, t - t') tlj{k, t') = -r]{k, t), [5.1] 



Tr{k, t-t') = - [ dujlmt{uj, k)e-'^^'-''^ (5.2) 

In order to study the initial value problem as a perturbative expansion, it proves conve- 
nient to write this self-energy in the following form 

j:^{k,t-t') = ^F{t-t',k) 

Fit -t',k) = - f —lmt(uj, ^)e-'"(*-*') (5.3) 

TT J UJ 



Using the form of the self-energy given by eqn. (|3.17| ) we see that the function F can be 
written as 

F{t - t', k) = - [v+{k)F_{t - t', k) + V-{k)F+(t - t', k) 

F±{t -t',k) = - f —lmt±{u, ^)e-'"(*-*') (5.4) 



We now integrate by parts the non-local term in ( ^.Ij) and use the fact that the adiabatic 
switching-on of the external Grassman current leads to a vanishing time derivative for the 
fermionic expectation value for t < ||3^. Upon switching off the external current at t = 
the effective Dirac equation of motion for the induced expectation value for t > becomes 

i-fo— - 7 ■ fc ^{k, t) + F(0, k)^{k, t) - dt' F{t - t' , k) ^{k, t') = 0, (5.5) 
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A perturbative expansion of this equation is obtained by writing 

Fit, k) = g^F^^\t, k) + g^F^^\t, k) + ■ ■ ■ 

i;{k,t)= V^(°) {k,t) + (7=^V^(2) (^^ ^) + ^4^(4) ^k,t) + --- 



(5.6) 



replacing these expansions in ( ^.51 ) leads to a hierarchy of equations, the first two terms of 
which are given by 



'ilo^^--fk\tlj^'\k,t) = 



(5.7) 



27o|- - 7 ■ fc ) V^'^ {k, t) = -F(2) (0, k)i;^'^\k, t) + f dt' F^^) (t - t', k) ^j^^^ [k, t') (5.8) 



(5.9) 



The solution of the zeroth order equation (|5.7|) is given by 

^(o)(fc,t) = ^ [BfW{k)e-'''' + DfW\-k)e'^' 



while the second order equation ( p.8[ ) can be solved in terms of the retarded free field Green's 
function 



S{k,t- t') = — [p_(fc)e-*'=(*-*') + P+(^)e*'^(*-*')] e(t - t') 
by proposing the form 

^(2)(fc,t) = Y. [B'il{t)U\k)e-^'' + Df^{t)V^{-k)e^ 



ikt 



(5.10) 



(5.11) 



>{2), 



1*{2) 



with the coefficients Bj^ (t) , Dj^ (t) being slowly time dependent. We focus on the time 



k,s 



k,s 



evolution of an initial state of particles near the Fermi surface, i.e, we set D 
k ^ kp- A straightforward computation leads to the following result 



*{0) 
k,s 



and 



B'i'^\t)=zBf 



k,s 



t— / — Im T,^(uj,k) 

TT J UJ 



k [ duj ^ ~ , -, 1 
+ _ / _ImS_(^,A;)- 

IT J UJ [UJ — k) 



sin[(a; — k)t] 



Bf'\t)=B^'^ 



k,s 



k,s 



(uj-k) 

■^/^ImE_(u;,A;) [^-^"^[("-/^^" 
n J UJ [uj — ky 



Bff{t)=iBftRet.{uj = k,k) 



k.s 



k.s 



(5.12) 
(5.13) 



Using the formulae in the appendix of ref . [^ , we find that in the limit t ^ oo the two 
terms in ( p. 121 ) can be combined to yield 



(5.14) 
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?(2,6) 



The contribution 5^ ' (t) in the t ^ oo and k ^ jj, can be understood by writing the 



integral in eqn. ( |5.13| ) as follows 
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duj Im S_(u;, k) 



cos[(u; — k)t] 



[00 
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UJ 



uj — fi , k = k — fi 



(5.15) 



Using the formulae found in the appendix of ref. p8[ we find for A; <^ M and Mt ^ 1 



that 



B 



' k.s 



A 



k,s 

ttA \k 

,2 



9 



247r2 



-Tkt-2\ln{Mt) 
Im S_(cj = k) 

M = Me^ 



and combining the above results leads to the final lowest order perturbative result 
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(0) 
k,s 



2Xkln 



i6kt - Tkt - 2A In(Mt) 

k 

—Re S_(u; = k) 



M 



(5.16) 
(5.17) 
(5.18) 

(5.19) 
(5.20) 



with Ti^{uj) given by ( [4.30| ) and 7 is the Euler-Mascheroni constant. We note that the 
damping rate Tk given by eqn. ( |5.17|) coincides with the result obtained in ref. |25|, but 
the logarithmic term is a novel result of our analysis and is the dominant contribution for 
particles at the Fermi surface. 

Then up to second order we find that the coefficient of the positive energy spinors in the 
expectation value of the spinor field is given by 



>(0) 



Bk s (t) = Bf 1 - t6kt - Tkt - 2A In(Mt) 



(5.21) 



The linear and logarithmic secular terms in the solution invalidate the perturbative ex- 
pansion at very long times. The dynamical renormalization group introduces a systematic 



resummation of these secular terms as described in ||38|. The implementation of the dy- 
namical renormalization group begins by recognizing that the terms in the square bracket 
in ( |5.21|) can be interpreted as a renormalization of the coefficient Sl . Thus following the 
procedure detailed in fSH] we introduce the renormalization constant Z{t) in the form 



B 



(0) 
k,s 



BkAr)Z{r) 



Z{t) = 1 + \zi{t) + 



(5.22) 



with r an arbitrary time scale. The coefficient zi (r) is chosen to cancel the secular terms at 
the time scale t = r. Choosing this scale so that the perturbative expansion is still valid we 
find that the improved solution is given by 



BjiM = Btsir) 



i5k{t 



r 



Tkit 



2Aln(- 



r 



(5.23) 
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since the scale r is arbitrary the invariance of the solution on the choice of this scale, i.e 

dBUt) 



dr 







leads to the dynamical renornialization group equation, which to order A is given by 



dr 



Btsir) 



2()fc + 1 fc H 

r 







(5.24) 



(5.25) 



Finally we choose the arbitrary scale r to coincide with t in the solution of ( ^.25|) ^ 
leading to the renornialization group improved time dependent amplitude 



BuM = BiMe 



-i5k{t-to)p-Vk{t-to) 



2X 



(5.26) 



Therefore the time evolution of the wave function of quasiparticle states near the Fermi 
surface is given by 



i^k^kAt) = ^fc.fc,(to)e-^('=+^^)(*-*°)e-r^(*-*") 



n2A 



(5.27) 



We see that the oscillation frequency k + Sk coincides with the quasiparticle dispersion 
relation ( 4.68| ) to lowest order in A while the "damping rate" F^ coincides with the quasi- 
particle width (|4.72| ) to lowest order in A. The novel aspect is the power law relaxation 
with anomalous dimension which cannot be extracted by computing the damping rate in 
an equilibrium formulation. Thus we are led to one of the important results of this arti- 
cle: quasiparticle excitations with the Fermi momentum relax with a power law with an 
anomalous dimension. 



B. Time evolution from the resummed spectral density 

1. k = 

We now use the renornialization group improved form of the propagator for particles 
near the fermi surface to obtain the real time evolution from the inverse Fourier transform 
of (|3.18|) . An initial value problem is obtained by introducing an adiabatically switched-on 
external Grassman source that vanishes at t = whose Fourier transform is given by 



r]{uj,k) 



^7° 



uj-iO+ 



E [B^k,sU%k) - Diy\-k) 



(5.28) 



An initial particle state is prepared by choosing D^ s~^ with V-{k)ri{u!, k) = 2j'^ri{u, k). 
The Fourier transform is obtained by writing e"*"^* = e-«A»*e-«'^* and performing the analytic 
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continuation ( [4.37| ). For k = fi the long time behavior is dominated by the region K ^ for 
which the analytically continued propagator is given by ([4.47|) , we find 

(5.29) 

The power law can be extracted by a simple scaling of the integration variable K = z/t 
and is clearly a consequence of the anomalous scaling behavior of the propagator. Thus 
the time evolution obtained directly from the Fourier transform of the (Euclidean) renor- 
malization group improved propagator confirms the result from the dynamical (real time) 
renormalization group for excitations with Fermi momentum {k = 0) and relates the power 
law fall-off to the anomalous scaling dimension. 



2. k^O 



The Fourier transform with the improved propagator ([4.66|) is more complicated because 
of the branch cut and the complex pole. Although obtaining an exact expression for the full 
integral is a complicated task, we can confirm the exponential relaxation in the narrow width 
approximation (weak coupling) from the contribution from the complex pole (|4.75| ). The real 
and imaginary parts of which are given by ( |4.68| ) and ( [4. 721 ) respectively in the narrow width 



approximation. Keeping only the contribution of the complex pole to the Fourier transform 
for A; 7^ we find 

^Pik, t) '''^' ^ik, o)e-*'^*e-''^^(^)*e-^(^)*. (5.30) 

The exponential relaxation confirms the result of the dynamical renormalization group 
since the narrow width approximation implies A << 1 for which the power law arising from 
the cut contribution is subleading. 

Thus the study of both cases /c = and A; 7^ via the Fourier transform of the renor- 
malization group improved propagator confirms the results obtained directly in real time 
via the dynamical renormalization group. As was emphasized above, the most important 
novel feature is the power law relaxation with anomalous dimension for quasiparticles with 
Fermi momentum as well as the anomalous scaling dimension for the quasiparticle dispersion 
relation, width and residue. 



VI. FROM QED TO QCD 

Although we have focused our calculations on the case of QED, it is straightforward to 
extrapolate the results to QCD to the same order in the HDL approximation. The reason 
is that the polarization tensor for gluons (photons) to this order is given by the quark loop. 
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In the cold non-abelian theory, the gluons give only the vacuum contribution while the 
quarks dominate the polarization tensor for large density. This is unlike the case at finite 
temperature where the gluon loops give a contribution of the same order as that of the quark 
loop. Furthermore, again unlike finite temperature, a magnetic gluon mass is not expected 
to arise because the gluons do not have any infrared divergence. The important changes 
that are required to extrapolate the results from QED to QCD are the following: 

a): The gluon mass scale (|4.27| ) that enters in the longitudinal and transverse spectral 
densities ([4.26|) and ( [4.25|) respectively now becomes 



QCD - 4^2 V 2 

with Np the number of fiavors of quarks in the fundamental representation. 

b) : The effective coupling in the self-energy of the quark now includes the trace over the 
color matrices resulting in that for QCD 

^QCD - 24^2 2iV, 

Thus the results obtained for QED to leading order in the HDL approximation can be 
extrapolated to QCD via the replacement M -^ Mqcd ', A — *> Xqcd- 

Vertex corrections: an important question that must be addressed is that of vertex 
corrections. In leading order in HDL, QCD and QED are similar because the quark loop gives 
the leading contribution to the vector boson polarization tensor, therefore the only relevant 
vertex is abelian. Since the vertex is related to the self-energy by the Ward identity one 
would expect logarithmic corrections arising from the vertex. However the results obtained 
from a detailed study in |^ show that this is not the case. The basic point of the argument 



is the following, the small momentum limit of the vertex, which is the relevant limit for long 
wavelength gluons, is related to the quark self-energy as 

n,«0) = fM (6.1) 

Since the leading behavior near the Fermi surface is 

E(c<j, k) (X {u — ji) In \uj — jj] (6.2) 

only the time component of the vertex would in principle lead to an infrared divergence, 
however, this vertex corresponds to the exchange of a longitudinal gluon which is Debye 
screened and does not lead to an infrared divergence. The spatial components of the vertex 
are infrared finite, and do not change the leading infrared divergence. 

The main conclusion extracted from this argument and based on the detailed analysis 
of reference ||2^ (to which we refer the reader) is that the infrared divergences in leading 
order in HDL do not receive vertex corrections. Hence we expect that the resummation of 
the leading infrared divergences provided by the renormalization group either Euclidean or 
in real time do not receive vertex corrections. 
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VII. CONCLUSIONS, AND CONJECTURES 

Our goal in this article is a systematic study of the breakdown of the Fermi liquid de- 
scription of the normal state of cold and dense QED and QCD. As was recognized before, 
the exchange of magnetic gluons and photons that are only dynamically screened by Landau 
damping lead to infrared divergences in the quark self-energy for particles near the Fermi 
surface. The main aspect of this article is the recognition that these divergences are akin to 
those arising in a critical theory near its upper critical dimensionality. We then implemented 
both the Euclidean and the dynamical renormalization group to provide a non-perturbative 
resummation of the leading infrared singularities to analyze equilibrium and non-equilibrium 
aspects of the single (quasi) particle states near the Fermi surface. The non-perturbative 
resummation of the single particle propagator leads to a novel description of the spectrum of 
quasiparticle excitations with momentum near the Fermi momentum summarized as follows 

(i) The quasiparticle pole, width and residue (wave function renormalization) are given 

by 



Up{k) = sign{k) |A;|M-^^ cos[7rA] 



r(A;) = Zp[A;]|A;|sin[7rA] 

M 



Zp[k] 



(1-2A) 
oj = ui — kp ; k = k — kp 



respectively. Therefore the residue of the "quasiparticle pole" and the "quasiparticle 
width" vanish near the fermi surface as 



Zp[k] oc |A; — kpl ^ 
rCk) oc \k-kF\~ 



2\ 
-2A 



(7.1) 
(7.2) 



and the group velocity of quasiparticles with Fermi momentum vanishes indicating a 
collective backflow. 

The real time evolution of the single quasiparticle wave function for states near the 
Fermi momentum is given by 



ipk^kp(t) ^^fc«fcf(to)e' 



-J[fcF+cip(A;)](t-to)g-r(fc){t-t(,) 



n2A 



(7.3) 



Therefore quasiparticles with the Fermi momentum relax with a power law determined 
by the anomalous scaling dimension revealing that the physics near the Fermi surface 
is similar to that of a critical theory. 
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(ii) The single quasiparticle distribution function is continuous near the Fermi momentum, 
i.e. the "jump" discontinuity in the Fermi-Dirac distribution vanishes as a consequence 
of the vanishing of the quasiparticle residue at the Fermi momentum. We find 

sinivrAl (k — kp) ,^,, , so ._ ,s 

The vanishing of the discontinuity of the distribution function at the Fermi surface is 
the hallmark of the breakdown of Fermi liquid theory. 

There are some remarkable similarities between these features of cold dense QCD and 
QED and those of quasi one dimensional Fermi systems with marginal interactions that lead 
to a description as Luttinger liquids [^,^]. These systems also feature spectral densities 



and correlation functions with anomalous dimensions which are non-universal and depend 
on the couplings p9| , ^ . The fermion propagator that results from one gluon exchange 
(dynamically screened) is also similar to that conjectured to describe the normal phase of 
high temperature superconductors in terms of marginal Fermi liquids |^,J9|. Even if these 



similarities are just coincidental, it is possible that a body of results in the condensed matter 
literature could be useful to understand the fundamental aspects of the normal phase of 
QCD at high density. 

Potential impact of the results: As we mentioned in the introduction, the properties 
of the normal state not only determine the thermodynamics and transport properties but also 
influence those of the superconducting state, since superconductivity appears as an instability 
of the normal state towards pairing interactions. Therefore the results obtained in this article 
that pertain solely to the normal phase could bear on characteristics of the superconducting 
state as well as transport in the normal (and possible in the superconducting) phase. 



i): In reference |27] it was found that the same type of infrared divergences that lead to the 
breakdown of Fermi liquid theory are responsible for a substantial decrease of the su- 
perconducting gap. In the calculation of reference |2^, these divergences are manifest 
through the wave function renormalization in the equation for the gap. Furthermore, 
recently lifetime effects (of the normal state quasiparticles) were shown to lead to a 
further decrease of the superconducting gap [0. In this article we have shown how the 
wave-function renormalization logarithmic divergence sums up to produce anomalous 
dimensions, which change the quasiparticle description. The spectral density features 
anomalous dimensions for frequency and momenta near the putative Fermi surface as 
well as important modifications in the quasiparticle dispersion relations, lifetimes and 
residues. It is therefore a relevant question to assess the potential corrections of the 
resummed spectral densities for the normal quark propagators to the superconduct- 
ing gap. This entails solving the gap equation or alternatively the Dyson- Schwinger 
equation but with the normal quark propagators replaced by those obtained from the 
resummed spectral functions. We are currently studying this possibility and expect to 
report on our findings in the near future. 
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ii): As mentioned in the introduction the most relevant physical setting for cold-dense 
QCD is that of astrophysical compact objects, in particular protoneutron stars from 
Type II supernovae collapse or neutron stars (or pulsars). An important observational 
aspect that could yield information from the core of neutron stars is cooling, which is 
studied through (soft) X-ray emmision |]SD|. The cooling equation 



dE dT 

relates the neutrino and photon luminosities to the specific heat C. For normal quark 
matter the neutrino emissivity is dominated by the direct quark Urea processes d — > 



u eV ] u e ^ du [^. A norma/ weakly interacting degenerate Fermi gas has a typical 
specific heat linear in temperature C{T) = CqT where the coefficient Co depends on 
the Fermi liquid properties. In particular the wave function renormalization (residue 
at the quasiparticle pole) renormalizes the free field (free Fermi gas) value of Cq. Since 
the quasiparticle residue vanishes, signaling the breakdown of Fermi liquid theory, we 
conjecture that the linear law will be replaced in a manner similar to that of Kondo- 
type systems [Q C oc ln(T). If a superconducting instability introduces a gap for all 
quarks, then the specific heat will have a typical behavior Ci e~^*^"^^/^ (for T < T^) 
where A(T) is the superconducting gap and Ci depends on the density of quasiparticle 
states near the Fermi surface which again receives non-Fermi liquid renormalization 
corrections. If the neutrino and photon emissivity of quarks is suppressed because of 
the presence of a superconducting gap, then the electronic specific heat will become 
relevant, in which case again the non- Fermi liquid corrections to the specific heat will 
be important since there are no superconducting pairing instabilities in QED. If some 
of the quarks are un-gapped these are normal and will receive renormalization of their 
Fermi liquid behavior. 

Understanding the potential corrections of the breakdown of Fermi liquid theory both 
for quarks and leptons could therefore lead to a deeper understanding of cooling in neutron 
stars and therefore justifies further studies along these lines. 

Our analysis was based on a perturbative approach, therefore its range of validity is re- 
stricted to asymptotically large densities so that the effective coupling as(/^)/6vr ^ 1. Using 
the running of the QCD coupling constant, the weak coupling condition implies chemical 
potentials and baryon densities orders of magnitude larger than those available at the core of 
neutron stars. Thus our studies for the properties of the normal phase, along with those of 
color superconductivity must be taken as indicatives of novel and potentially relevant phe- 
nomena, but must be carefully extrapolated to the realm of baryochemical potential (baryon 
density) relevant to quark matter in neutron stars, requiring alternative non-perturbative 
techniques. 
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APPENDIX: REAL-TIME PROPAGATORS 

1. Dirac Fields 

In this appendix we summarize the various real-time propagators used in this article. 
The fermion propagators are defined by 

(M/«(f , t)^\x', t')) = 1 1 -0^,sfit, t') e^-^-^\ 

S^^t,t') = S^{t,t')9it - f) + SS{t,t')9it' - t), 
S^-it,t') = S^{t,t')9it' -t) + S^it,t')9{t - t'), 
Sf{t,t') = Sl{t,t% (Al) 

where a, b = ±. 

In an equilibrium situation the propagators can be written in terms of spectral densities 
as follows 

tS^{t,t') = ldqoP>{qo,k)e-'''^'-''^ 

p>(go, k) = p{qo, k){l - Nf{qo, k)) 
p<{qo,k)=p{qo,k)Nf{qo,k) (A2) 

Where the Fermi-Dirac distribution functions for particles and antiparticles, for the case 
under consideration of finite chemical potential p and zero temperature are given by 

Nf{q,) = e{p-qo) (A3) 

Nf{qo) = e{-p-qo) (A4) 

For free fields the fermion Wightman functions are given by 

S|(t, t') = -^{(^ + m)[l- NAuj^)]e—^^'-''^ 
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+ 7o(^ - m)7oiVp(cc;g) e*-s(*-*')}, 

+ 7o(y^ - m)7o[l - Nf{u;^)Y^^^^'-''^], (A5) 



with i^ = (u;^, /c), /^ = 7°C(;^ — 7 ■ A; and c^;^ = \J k'^ + rn?. For massless fermions (the case 
under consideration) 



'^■ + '">] =T>_(j)=7°-7-* 



fe / m=0 



with the properties 



^k ) m=0 



(A6) 



V.{k)V+{k) = 2-i'^V+{k) 

V+{k)V- {k) = 27°P_ (k) (A7) 



2. Gauge Fields 

a. Spatial components 
The transverse photon propagators are defined by 



,bf;^ i.>\\ _ ■ I Ji_S_n"P (+ +'\ -vi^ ( Fr\ ^i<i<^-^') 



Gt.i^^t') = GrAt.t'nt - t') + ^^,,(t,t')^(t' - t), 

^--(t,t') = gi^{t,t')e{t' -t) + gi,{t,t')e{t - f), 

^±;(t,t')=^|,,(t,t'), (A8) 

where Vq^{q) = 6''^ — g^g-' is the transverse projector and the photon Wightman functions 
can be written in terms of a spectral representation as follows: 

OT,git,t') =tJdqoPT{qo,q)nB{qo)e'"'''^'~''\ (A9) 

where priqo, q) is the spectral density and nB^qo) is the Bose-Einstein distribution function. 
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At zero temperature 

riBiqo) = -©(-go) ; 1 + nsigo) = 0(go) 
For free fields at zero temperature we find 



^TyM') = ^e- 



-iq{t-t') 



QUt^t') 



J_ iq(t-t') 

2g 



In the HDL approximation, the spectral density priQo, o) is given by 
PT(go, q) = sgn(go) Zxiq) 6[ql - u^iq)] + Priqo, q) 0{q^ 



ql) 



47r2 fc V k^J 



{ 



uj'^ — k"^ 



47r2 



2a^ 
fc2 



+ i(i-#)iHf^I1} +['^i(i-#) 



2- 



and ^tI?) is the residue at the pole of the collective excitation |^ . 



(AlO) 



(All) 



(A12) 



h. Time component of the photon propagator 

The time-time component of the photon propagator describes the instantaneous Coulomb 
interaction including screening corrections and are given by 

ht^i ,>^^_, f d'q gab^^^^,yqiS-S')^ 



{Al{x,t)A',{x\t'))=i 



(2vr) 



^++(t,t') = -5{t - 1') + Gl,{t,t')e{t - 1') + gut,t')e{t' - 1), 
^^,-(t,t') = -\6{t - 1') + gi^{t,t')e{t' -t) + gut,t')e{t - f), 

with the Wightman functions expressed in terms of the spectral density pi as 
Glqit, t') = -^j dqo hiqo, q) [1 + ns^qo)] e-^''^^*"*'), 
Qlqit. t') = -ij dqo pLiqo, q) UB^qo) e'^'^^^'-'\ 
For free fields piiqo, ?) = and 

Qt:,{t,t') = \6{t-t'), 

:y qZ 

,y qz 

gSit,t') = o 
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(A13) 



(AM) 



(A15) 



In the Hard Dense Loop (HDL) approximation the spectral density Pl{%i q) is given by 
Pl(?o, q) = sgn(go) ^L(g) S[ql - ujUq)] + f^iigo, q) 0{q^ 
h{qo,q) -- 



ql) 



27r^ q 



,2 ^^(2 -20 In 



q' + 



g+go 
q-Qo 



I rsfiiigois 



(A16) 



where ^^(g) is the plasmon (longitudinal photon) pole and Zi^q) is the corresponding 
residue 
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FIG. 1. p_(^,^) for A = 0.1; ^ = 0.2 
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FIG. 2. p_(^,f)forA = 0.01;A = 0.2 
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